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ImpZ: a new photometric redshift code for galaxies and 
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ABSTRACT 

We present a combined galaxy-quasar approach to template-fitting photometric red- 
shift techniques and show the method to be a powerful one. 

The code, ImpZ, is presented, developed and applied to two spectroscopic redshift 
catalogues, namely the Isaac Newton Telescope Wide Angle Survey ELAIS Nl and 
N2 fields and the Chandra Deep Field North. In particular, optical size information 
is used to improve the redshift determination. The success of the code is shown to be 
very good with Az/(1 -I- z) constrained to within 0.1 for 92 per cent of the galaxies in 
our sample. 

The extension of template-fitting to quasars is found to be reasonable with Az/(1-|- z) 
constrained to within 0.25 for 68 per cent of the quasars in our sample. Various tem- 
plate extensions into the far-UV are also tested. 

Key words: galaxies: evolution - galaxies :photometry - quasars: general - cosmology: 
observations 



1 INTRODUCTION 

Photometric redshifts are a powerful statistical tool for stud- 
ies of the evolutionary properties of galaxies, in particular 
of faint galaxies, for which spectroscopic data is hard or im- 
possible to obtain. 

Photometric redshifts are faster to measure than spec- 
troscopic redshifts and can be applied to much fainter 
magnitudes since the bin sizes are larger («1000A vs. 1- 
2A). There is how ever a trade-off with redshift precision - 
iHogg et al.l 11 19981) found that photometric redshifts can be 
predicted with an accuracy of Az ~0.1 (0.3) in z for 66 per 
cent (99 per cent) of the sources examined. 

The bulk of the photometric redshift identification is 
carried out using the broadband continuum shape and pres- 
ence and/or absence of spectral breaks (like the 4000A 
Balmer Break) or the onset of the Lya forest and Lyman 
limit which enter opti cal waveband s at high redshift (the 
Lya forest effect, e.g. iMadad llQAiil and ISteidel fc Sareend 
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119871) . Although methods based on training sets such 
as polynomial fitting (e.g. IWang. Bahcall_fc_Tu rner 199^ 
and artificial neural networks (e.g iBall et aL 2004 and 
iTagliaferri et al]|2002^ have had some success in determin- 
ing photometric redshifts, we decide to utilise the template- 
fitting procedure (but with an element of training in that 
the templates have been adapted to improve fits). This is 
because there is a relative paucity of galaxies with both spec- 
troscopic redshifts and photometry in the bandpasses used 
in the catalogues under study here, particularly at higher 
redshifts, making the construction of a training set unfea- 
sible. Empirical methods are also hard to apply outside of 
the boundaries in which they were defined - such as the 
redshift distribution of the training set or the photometric 
bands use. 



In the template-fitting method fe.g. lSawicki. Lin fc Yed 

I1997I : .Giallongo et al. 1998) the observed galaxy fluxes, 
fi^" in the i"" band, are compared to a library of reference 
fluxes, f*'^™-P (z,T), where T is a set of parameters that 
account for the template galaxy's morphological type, age, 
metallicity and dust. We then fit the observed fluxes to the 
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library fluxes using x^ minimisation. As well as deriving 
redshifts, the procedure produces information on spectral 
(template) type, although this is less robust due to the 
degeneracies in the parameter space (see H6. 3.111 . A way of 
brea king these de generacies is to use Bayesian probability 
(e.g. l,Tavnesll200ril to weight the solutions based on a prior 
knowledge of the expected population distributions. Appli- 
cation of Bayesian methods to photome tric redshift s has 
been pr e sented in, for example, Kodama et alJ (|l99{y and 
iBenitea (|2000|). Usually such applications use priors such 
as the expected redshift distribution of the sample but this 
naturally suppresses unbiased information on the true red- 
shift distribution. Here Bayesian methods are implemented 
using absolute magnitude limits and extinction distributions 
which it is hoped retain the power of Bayesian methods 
without unduly influencing the underlying science (see t|2.2ll . 

One important refinement is the consideration of 
extinction for galaxies. iMadaul (|l999|) used a single dust ab- 
sorption correction of Ai5oo = 1.2 mag, except for galaxies in 
the redshift range 0.75-1.75 where the equivalent extinction 
at 2800 A was used. Galaxy evolution mo dels such as those 
of iLe Borene fc Rocca-Volmeranea (|2002h include evolution 
of dust extinction with time. In order to allow for variation 
in extinction from galaxy to galaxy, extinction needs to 
be solved as an ad ditional free parameter to redshift. In 
ISteidel et alJ (Il999h dust absorption was corrected for by 
assuming that colour deviations in their sample galaxies 
were entirely due to dust absorption, based on the empirical 
relation between far-infrared emissi on and the observe d 
UV spectrum s lope, as derived bv [ Meurer e t al. fl9 9J|). 
More recently, Ixhompson. Wevmann fc Storric-Lombardil 
i200ll) used Spectral Energy Distribution (SED) template- 
fitting photometric redshift techniques in the deep 
NICMOS northern HDF, fitting extinction as a pa- 
ra meter ranging from E{B — V)=Q to 1.0. The study 
of iBolzoneUa. MiraUes fc Pelldl ll200(i) found that the 
inclusion of A„ as a free parameter caused significant 
increases in aliasing. In a similar technique developed in 
[Rowan- Robinson (2003a), hereafter RR03, these aliasing 
problems were reduced by setting several A,, priors (see 

am- 

In this paper the SED template-fitting set out in 
RR03 is refined and applied to two spectroscopic redshift 
galaxy sampl es from the Europ ean Large-Area ISO Sur- 
vey (ELAIS: loiiver et alJlJOOol) Nl and N2 fields of the 
Isaac Newton Telesco pe Wide-Angle Survey (INT WAS; 
iMcMahon et al.ll200ir i. a part of the INT Wide-Field Sur- 
vey (INT WES), and also the Chandra Deep Field North, 
CDEN) for validation purposes. The eff'ect of non-zero A„, 
different SED templates and various template extensions 
into the far-UV are explored, as is the inclusion of several 
different quasar templates and the applicability of template- 
fitting techniques to quasar-like sources. In JUthe photomet- 
ric redshift technique is set out. In i|21the various templates 
and extensions to the UV are discussed. The ImpZ code is 
then applied to two spectroscopic redshift galaxy samples 
in [21 The results of this validation are given in JH Error 
analysis is discussed in ^\ Discussions and conclusions are 
presented in ^and ^ 

The application of the ImpZ code to the entire re- 



calibrated ELAIS Nl and N2 fields from the INT WAS and 
investigations into the evolution of extinction and star for- 
mation rates (SFR) will be presented in a companion to this 
paper; Babbedge et al. (2004; in prep.). 
Note that for these investigations the fiat, J1a=0.7 cosmo- 
logical model with Ho=65km s~^Mpc~^ is used. 



2 METHOD 



2.1 



X 



^ analysis 



The template-fitting procedure is as follows: the observed 
galaxy magnitudes are converted for each i*'' photometric 
band into an apparent flux, /"'"'. Equivalently, this recon- 
structs the target galaxy's SED at a very low spectral resolu- 
tion by sampling the luminosity at the effective wavelength 
of each photometric band. The observed fluxes can then be 
compared to the template fluxes, f^''"^^ (z,T), computing the 
reduced x^ , xled as 



2 
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(1) 



where N is the number of photometric bands, D is the num- 
ber of degrees of freedom, ai is the observational uncertainty 
in the i*'' band (hence the solution is weighted by the flux 
errors, however to avoid excessively high weighting by very 
high signal-to-noise observations, a minimum flux error is 
set, typically 0.5 per cent) and s is a normalisation factor to 
minimize x^ for each template. 
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(2) 



Note that if there is a detection in just one band then fitting 
is not attempted. 



2.2 The ImpZ code 

This code builds on the technique presented in RR03, ex- 
tending it to include the correct treatment of CCD response, 
filter transmission characteristics and the statistical effect of 
Inter-Galactic Medium (IGM) absorption (Ly g, Ly/3, Ly7, 
Lyi5, and Lyman continuum), as set out in iMadau et alJ 
|l99^ (see Fig. 0. Correct treatment of the effects of the 
IGM is needed because it is possible to mis-interpret the 
Lya forest effect as the intrinsic Lyman break and because, 
particularly at high redshifts, it will imprint its own recog- 
nisable feature onto the SED. Although the treatment of 
IGM absorption is only based on the average accumulated 
absorption, the rms fluctuations away from the mean can 
be expected to be small once integrated though a broad 
band pass (Press. Rybicki fc Schneid er 1993). At high red- 
shift, |MassaxottijJovirw_&_^uzzoni| (1200111 have shown that 
correct treatment of internal dust reddening (the Interstel- 
lar Medium) and IGM attenuation are the main factors in 
photometric redshift success. The effect of internal dust red- 
dening for each galaxy is already incorporated in the tem- 
plates (Table |5| a nd alterable via fltting for A„, using the 
reddening curve of lSavaee fc Mathij (Il979l) . Observed fluxes 
are compared to template fluxes for 0.01<logio(l-l-z)^0.90, 
equivalent to 0.02<z<6.94. 
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Table X. Common spectral lines and breaks in SEDs, and the redshifts where they enter and leave the filters used in this study, where transition is defined 
by the points w^here filter transmission drops below 50 per cent of the maximum filter transmission. Redshifts greater than seven are not considered. 



INT WAS filters 




Lyman 


Limit {912A) 


Lyo 


(1216A) 


Balmer Break (4000A) 




doublet (4983A) 


Ho 1 




Filter 


0[///] 


[6563A) 




enters 


leaves 


enters 


leaves 


enters 


leaves 


enters 


leaves 


enters 


leaves 


U 


2.5 


3.3 


1.6 


2.2 


- 


- 


- 


- 


- 


- 


9' 


3.5 


4.9 


2.4 


3.4 


0.03 


0.35 





0.08 


- 


- 


r' 


5.0 


6.6 


3.5 


4.7 


0.38 


0.73 


0.10 


0.38 


005 


0.07 


i' 


6.7 


NA 


4.8 


6.0 


0.75 


1.1 


0.40 


0.71 


0.07 


0.30 


Z 


NA 


NA 


6.0 


6.9 


1.1 


1.4 


0.71 


0.93 


0.30 


0.46 










Barger 


et al. (2000) filters 












Lyman 


Limit (912A) 


Lya 


(1216A) 


Balmer Break (4000A) 




doublet (4983A) 


Ha 1 




Filter 


OlIII] 


;6563A) 




enters 


leaves 


enters 


leaves 


enters 


leaves 


enters 


leaves 


enters 


leaves 


B 


3.3 


4.5 


2.2 


3.1 





0.25 





0.003 


- 


- 


V 


4.5 


5.6 


3.1 


3.9 


0.25 


0.50 


0003 


0.20 


- 


- 


R 


5.6 


6.7 


3.9 


4.8 


0.50 


0.75 


0.20 


0.40 





0.07 


I 


6.9 


NA 


4.9 


6.1 


0.80 


1.2 


0.44 


0.73 


0.10 


0.31 


Z 


NA 


NA 


5.5 


NA 


1.1 


1.5 


0.69 


1.0 


0.28 


0.52 


HK', (H) 


NA 


NA 


NA 


NA 


2.7 


3.5 


2.0 


2.6 


1.2 


1.7 


HK', (K') 


NA 


NA 


NA 


NA 


3.9 


4.8 


2.9 


3.6 


2.0 


2.5 



The following parameters and cuts were used: 

(i) in order to reject unphysical fits only tfiose that give 
absolute B-band AB magnitudes (Mb , 4400A) in the range 
[—22.5 — 2logio{l + z)] <Mb < —13.5 are considered for the 
6 galaxy templates. Having an upper envelope dependent on 
redshift was found to be the best way to suppress luminous 
outliers at low redshift whilst allowing more luminous galax- 
ies at higher redshifts, and this luminosity-redshift depen- 
dence agrees with the natural consequence of known strong 
luminosity evolution for galaxies (e.g. Lin et aL 1 9991) . For 
the AGN templates a range from —27 <Ms < —17.5 was al- 
lowed. These restrictions cut out excessive numbers of aliases 
at the minimum and maximum redshift. A number of differ- 
ent limits have been investigated but these were found to be 
the most effective. Originally, limits were applied to the / 
band but this failed to constrain the UV luminosity. Shifting 
the limits to the B band allows both the young and old star 
components of the SED to be constrained for the redshifts 
considered. 

(ii) Sources are defined as stellar or non-stellar, where 
this is a purely morphological property differentiating be- 
tween point-like and more extended sources. Those defined 
as stellar will have Active Galactic Nuclei (AGN) templates 
considered in addition to the galaxy templates, whereas non- 
stellar sources will only be fit by galaxy templates. This 
reduces the increased degeneracy introduced by including 
AGN templates. The procedure for splitting sources into 
stellar /non-stellar sub-groups is different for the INT WAS 
ELAIS spectroscopic sample and the CDFN sample, and is 
described for each in t|4.1l and H4. 21 respectively. 

(iii) sources that were saturated were removed since their 
photometry is poor. This reduces Zphot outliers. 

(iv) A prior expectation that a stellar source is more 
likely to be a QSO is introduced by minimising Xred + '^F{T) 
rather than x^ (a=4 here) for stellar sources. F{T) is a delta 
function such that F{T) = if the template, T, is an AGN 
template, and F{T) = 1 if the template is a galaxy tem- 



plate. Essentially it is a prior to prefer AGN solutions due 
to the morphology information contained in the class fiags 
- a weak Bayesian (e.g. .Tavnes 200(1) formulation - and was 
reached based on the results in >l5.1.4l and tl5.2.3l 

For the 'free' Ai, fitting option the following restrictions 
were used: 

(v) for the elliptical and AGN templates, At, can take the 
value only. The reason that A„ is set to zero for AGN is 
that since AGN are essentially a power-law, the additional 
inclusion of A„ gives too much freedom to the shape of the 
AGN template, and it was found that resulting degeneracies 
reduce the effectiveness of the photometric redshift tech- 
nique. 

(vi) for other templates, A„ can take the values 0.0 to 
1.0, in steps of 0.1. The maximum At, was chosen to be 
approximately t wice that of t he ty pical At, of galaxies at 
z ^ 4 found by ISteidel et al.l ^1999^ who derived a typical 
E{B — V) of 0.15. Note that the templates already include 
some extinction (see Table |2l so that the A„ of the solution 
is technically the difference between the actual value and 
the template's inherent value. 

(vii) no solution for A„ is sought if the reduced x'^, Xred) 
of the A„=0 solution is <1, or if there are less than 4 bands. 

(viii) a prior expectation that the probability of a given 
value of At, declines as | At, | moves away from is introduced 
by minimising Xred + P-^v rather than x^ (/3=2 here). This 
can again be viewed as a weak implementation of Bayesian 
methods. 



3 TEMPLATE SEDS 

The choice of how many templates to use in fitting is a 
crucial one. The choice of too many leaves the code with 
too much freedom, leading to large numbers of aliases and 
degeneracy. Similarly, too few and the code will be un- 
able to find accurate redshifts for real objects, something 
iBolzonella et alJ ||200C|) termed 'catastrophic' failures. 
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5000 



I 0,000 I 5,090 

Wavelength, A 



20,000 



25,000 



Figure 1. IGM transmission effect at rcdsliift 1, 2, 3, 4 and 5 on the observed ffux from a source. Below are sliown the INT WAS filters {solid lines), 
combined with the CCD response of the WFS camera, and also the filters (solid lines) used in the Barger et al. (2002) CDFN catalogue, combined vi'ith the 
relevant CCD (Suprime— Cam and QUIRC) responses. The additional effect of including IGM absorption is illustrated for the redshift 3 case by over-plotting 
the filter response when combined with the redshift 3 IGM transmission curve. It can be seen to effect U and g in the INT WAS filter set, and the B filter 
in the CDFN filter set (over-plotted as dashed lines). 



As in iMobasher et alj (Il99d) and RR03 the code uses 
six galaxy templates; E, Sab, Sbc, Scd, Sdm and star- 
burst galaxies. These six templates (or similar versions) 
have been found to provide a good low-rosolution repre- 
sentation of observed galaxy SEDs (Mobasher ct aOllQQa . 
RR03), and are used in preference to a large array of evolv- 
ing galaxy tem plates, as generated by ev olutionary codes 
such Eis that of iBruzual fc Chariot! lll993i) . since too much 
freedom is then available in fitting. The empirical templates 
used in RR03 were based on observations (starburst tem- 
plate adapted from ICalzetti fc Kinnevlll992r) and on obser- 
vations and colour synth esis (the remainders adapted from 
[Yoshii & Takaharal ll988r . In a method reminiscent of em- 
pirical techniques, these templates were adapted to improve 
photometric redshift results. 

The templates used in this work were generated by re- 
producing these original templates via spectrophotometric 
synthesis, in order to strengthen their physical basis. The 
original templates were convolved with filter transmission 
curves in order to create virtual datapoints. These data- 
points were then fit using a code which combines a given 
number of Simple Stellar Populations (SSPs), each weighted 
by a different SFR and extinguished by a different amount 
of dust. The procedure is based on the synthesis code of 
IPoggianti. Bressan fc Franceschinil J20ni| l and ha s been pre- 
viously applied to spectra in lBerta et alJ l|2003a|). The code 
minimizes the x^ obtained from comparing the new template 
to the datapoints. Minimization is based on the Adaptive 
Simulated Annealing algorithm. Details on this algorithm 
and on the fitting technique are given in ?. The templates 
are plotted in Fig. |21 and have been extended further into 
the UV regime as set out in 



3.1 SSP Populations 

The spectra of the SSPs have been computed with a Salpeter 
Initial Mass Functio n (IMF) between 0.15 and 120 solar 
masses, adopting the IPicklea l|l998f) spectral atlas and ex- 
tending its atmospheres ou tside its original range of wave- 
lengths with lKuruc3 Jl993h models from lOOOA to 50,000A, 



as described in iBressan. Granato &: Silva (|l993). Nebular 
emission is added by means of case B HII region mod- 
els computed through the ionization code CLOUDY of 
iFerland et alJ lll998fl . The adopted metallicity is solar. 

When generating the new templates, each SSP is 
weighted by a different SFR. For each SSP a uniform screen 
attenuation is adopted, using the standard extinction law 
(Rv=A.„/(£;(B - l/)=3.1. ICardeUi. Clavton fc Mathij|l989l) 
and adopting a different E[B — V). 

With 10 SSPs, we have a total of 20 free parameters 
(SFR and E{B - V) for each SSP), but the code automati- 
cally discards those populations which contribute less than 
one per cent to the total spectrum, at each wavelength. As 
a result each template is constructed from only a few SSPs. 
The E{B — V) values for the oldest populations (9th and 
10th) are constrained to be less than 0.2 (the characteris- 
tic extinction of the older quiescent stellar population in a 
sample of nearby galaxies is E{B — V)~ 0.1; ?). The to- 
tal E{B — V) for each template is obtained by comparing 
the non-extincted final spectrum {Sn) to the extincted ver- 
sion (Se) since in the V band the following relation holds: 
Se = 10"°* ^ Sn. For details of the populations and their 
contribution to each template see Table |5| 

Two fits were generated for the elliptical template — 
El and E2. El fits the small UV bump present in ellipticals 
from approximately 1000-2000A, a feature that is due t o 
emission from planetary nebulae JYoshii fc T akahara" l988l) . 
In order to prevent this UV bump being fit by young stellar 
populations, the SFR of the 3 youngest SSPs was set to zero. 
E2 consists only of the two oldest SSPs and fails to fit this 
UV bump. For this reason, E2 was not used as a template, 
but is plotted in Fig. |5|for interest. 

The template-extension into the UV is described in 13. 31 



3.2 AGN Templates 

As well as galaxy templates, the inclusion of a number of 
different AGN templates has been investigated (see Fig. |^ 
to allow the ImpZ code to identify quasar-type objects as 
well as normal galaxies. This is of particular interest for 
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Wavelength, A 



Wavelength, A 
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Figure 2. The 6 galaxy templates used. Dashed lines show the original RR03 templates {offset for clarity), solid lines shows the SSP generated versions, 
along with extension into the Far-UV (sub— lOOOA) as set out in section|3^ For the elliptical template two SSP generated fits are shown, which diverge below 
around 2000A. Line El. fits the UV bump that is due to planetary nebulae, whereas Line E2. does not. This is discussed in section |^^^J 



Table 2. The six galaxy templates and the SSPs that were used to create them. SFR is given in Mqi/t and normalised so that each template has a total 





SSP 










SED Template 














# 


Age (yr) 


At (yr) 


Elliptical(El) 




Sab 




Sbc 




Scd 




Sdm 




Sbrst 




SFR E(B-V) 


SFR 


E(B-V) 


SFR 


E(B-V) 


SFR 


E(B-V) 


SFR 


E(B-V) 


SFR 


E(B-V) 


1 


10*^ 


2 10*5 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


2 


3 10"^ 


4 10'' 


- 


0.100 


0.04 


- 


- 


- 


- 


4.097 


0.001 


41.17 


0.9465 


3 


810'' 


310" 


- 


- 


- 


- 


- 


- 


- 


- 


- 


35.67 


0.024 


4 


10^ 


4 10*' 


0.099 


- 


- 


1.461 


0.18 


- 


- 


- 


- 


- 


- 


5 


5 10'' 


6.2 10'' 


- 


0.332 


0.015 


- 


- 


2.576 


0.03475 


13.54 


0.004 


- 


- 


6 


10* 


1.25 10* 


- 


- 


- 


- 


- 


- 


- 


13.62 


4.685 


- 


- 


7 


310* 


2 10* 


- 


0.484 


0.12 


1.647 


- 


1.692 


0.12 


9.995 


- 


- 


- 


8 


5 10* 


3.5 10* 


- 


- 


- 


0.896 


0.003 


2.80 


0.032 


11.58 


0.0616 


- 


- 


9 


10^ 


1.25 10" 


0.396 


1.074 


0.017 


1.509 


0.137 


0.771 


0.1838 


1.111 


0.0818 


7.78 


0.0015 


10 


1.2 1010 


lOl" 


0.950 


0.854 


0.134 


0.747 


0.186 


0.756 


0.138 


- 


- 


- 


- 


A„ 


(young - 1st 3 SSPs) 


0.0 


0.12205 




0.0 




0.0 


0.0031 




1.38 


A„ 


(total) 




0.0 


0.2344 


0.3308 




0.271 


0.4303 




0.736 



application of ImpZ to the entire ELAIS Nl and N2 fields 
of the INT WAS in Babbedge et al. (2003; in pr ep.) since 
many ELAIS sources are expected to be AGN (Ol iver et alJ 
I2OOOI) . Fitting with AGN templates is only carried out for 
sources that have been defined as stellar, as described in 

am 

The last decade has seen a large rise in the number 
of optically selected high-redshift quasars and the existence 
of large sampl es of quasars (e.g the Sloan Digital Sky Sur- 
vey - SDDS - lYork et alJl200Cl) means that the derivation 



of photometric redshifts for quasars is gaining popularity 
as well as reliability. For galaxies the technique relies on 
the identification of continuum features such as the 4000A 
break - see Table Q for examples. For a featureless spec- 
tra a photometric redshift is far harder to determine, if at 
all. The majority of quasars can be characterised in the 
UV-optical region as a featureless continuum. Overlaid on 
this continuum however are a series of (mostly) broad emis- 
sion line featu res which contain a significa nt amount of fiux 
=a] ll99ll : iRichards et al.ll200J) . Also, at higher 
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redshifts the result of Lya forest absorption will imprint 
an additional redshift-dependent feature onto the contin- 
uum. Empirical redshift -colour relationships have been ap- 
plied to quasars, but such polynomial fitting is limited due 
to the nature of quasar colours, which can change rapidly 
across a small redshift range as an emission line passes in 
and out of a passband, or remain constant with redshift due 
to the featureless continuum - a relationship which is poorly 
reproduced by polynomial functions that vary slowly with 
redshift. An extension of this technique which implements a 
nearest neighbour (NN) estimator with the reference points 
derived f rom colours averaged over redshift bins was pre- 
sented in lRichards et al.l i2001f) . with around 70 per cent of 
predicted redshifts matching reasonably well with the spec- 
troscopic values (Az < 0.1 for 55 per cent and Az < 0.2 for 
70 per cent of their sample) . 

The alternative approach is a template-fitting 
procedure as used in several qu asar studies (e.g. 
iHatziminaoglou. Mathez fc Pella l200Cl) and adopted in 
this work. This approach requires a template or templates 
that can cover the observed range of spectral types, across 
a large enough wavelength span to be applicable for 
low and high redshift objects. Unlike galaxies, quasars 
have SEDs with similar power-law continua so the use 
of the mean spectrum of a sample of quasars is feasible. 
[Budavari et al. (2001) used the SDSS composite spectrum 
of Vanden Berk et al. (2001) . resulting in a slightly greater 
rms than tha t foun d in the empirical NN method of 
iRichards et alj 11200 ih . This suggests that the use of a 
single quasar template is not sufficient - perhaps one 
corresponding to broad a bsorption lines or Fell is also 
needed. A similar study bv lWu. Zhang fc Zhoul ll2004) with 
the SDDS composite achieved Az < 0.1 for 47 per cent and 
Az < 0.2 for 68 per cent of their sample. iBudavari et alJ 
(|2001j) went on to reconstruct 4 discrete templates in an 
iterative manner, gaining better results than with the 
empirical NN method. However, this is beyond the scope of 
this paper and is left as a future direction. 

In this study, several AGN templates were constructed. 
The optical-IR basis of these templates is outlined below, 
and extension to UV in 



(i) the SDSS median composite quasar spectrum 
iVanden Berk et alj 1200 J) covers a rest-wavelength range 
from 800 to 8555A and is constructed from 2204 quasars 
spanning 0.044^ z $;4.789. For wavelengths longward of 
Lya to optical wavelengths the continuum is well-fitted by a 
power law (f(A) oc A"-^) with a wavelength power-law index, 
Qa {ct\ = — (ctv-\-2)) , of —1.56. This is consistent with a num- 
ber of other works based on optical and/or radio-selected 
samples that have found power-law continuum indexes, as 
average values from spectra, from pho tometry, or using com- 
posite spectra, of —2 ^ oa ^ — 1 fe.g.lSch neider et alJl200lL 



i Brotherton et al.l EooH. ICarballo et al.l Il999l iNatali et al. 



Iigga . 'Zheng ct al; 'l99/. 'F rancij Il996l (Francis et alJll991 
and ICristiani fc Vio 1990) . The continuum blueward of the 
Lya emission line is heavily absorbed (the median redshift 
is 1.253) due to Lya forest absorption. Since this effect is 
strongly effected by redshift however and the median com- 
posite uses spectra across a broad range of redshift, little can 



be drawn from the absorption in that region of the SED. In- 
stead, the UV part of the SED will be treated separately 
- see t|3.3| The template has been extended to wavelengths 
longer than 8555A by uti lising the IR part of the average 
optical quasar spectrum of lRowan-RobinsonI (|l995|), slightly 
modified as set out in ?. This extends the template out to 
25^m with what is essentially a continuation of the contin- 
uum power-law. See the 'SDSS' line in Fig.|3 

(ii) in addition to the SDSS composite template, two sim- 
pler AGN templates were include d. These are based on th e 
mean optical quasar spectrum of iRowan-RobinsonI il995l) . 
spanning 400A to 25^m. For wavelengths longer than Lya 
the templates are essentially a\ w —1.5 power-laws, with 
slight variations included to take account of observed SEDs 
of ELAIS AGN (?). These two AGN templates are referred 
to as RRl and RR2, and are very similar to those used 
in ?. They differ from one another at wavelengths longer 
than 1/im where RR2 contains more flux. For wavelengths 
shorter than the Lyman limit, several UV behaviours were 
again considered - see H3.3I These templates can be seen as 
'RRl' and 'RR2' in Fig.^l 



3.2.1 A red quasar 

Should a much redder quasar be considered? There is a 
debate over the existence of a sign ifica nt population of 
red q u asar s (e.g. Webster et al.lll995l and IBrotherton et all 
l200l[) . iRich ards et al. (200^found that roughly 6 per 
cent of their homogeneously selected sample of 4576 SDSS 
quasars were red in comparison with even the reddest power- 
law continuum quasars and were probably dust reddened. 
In this work a red AGN template has not been included - 
some tests were done using the z ^ 2.216 FIRST J013435.7- 
093102 source from iGregg et alJ i2002f) . which is an ex- 
tremely dust-reddened lensed object with B— K^ 10 but 
it was found that it did not match any sources in our two 
spectroscopic redshift catalogues. The inclusion of a red 
AGN template is, however, expected to be powerful when 
optical data is combined with upcoming IR data from the 
SPITZER mission, particularly for the SP ITZER Wide- 
Area Infrared Extragalactic Survey (SWIRE jLonsdale et alJ 
l2003h which covers an area large enough (50 square degrees) 
to find a significant number of these rare objects. Appli- 
cation to SWIRE data will be considered in a future work 
(Babbedge 2004; PhD thesis, in prep.). 



3.3 Far UV treatment 

The Milky Way becomes virtually opaque at wavelengths 
between 100 - 912A due to absorption by neutral hydrogen. 
Hence it is difficult to be certain of the far-UV rest-frame 
emission of galaxies or quasars. Ground-based observations 
can only start revealing the sub-1200A regime for objects at 
z >2 due to atmospheric absorption of shorter wavelengths. 
Escaping our own atmosphere with space-based telescopes 
has allowed the far-UV and extreme-UV region to be ex- 
plored in more detail, but even then the action of the IGM 
and of neutral hydrogen in our own Galaxy makes the SED 
determination uncertain. In order to apply a template-fitting 
technique out to large redshifts it is important that the tem- 
plates extend to sufficiently short wavelengths. For galaxies 
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Figure 3. The various AGN templates that were investigated. Tlie solid 
line 'SDSS' is the mean SDSS quasar spectra (Vanden Berk ct al. 2001), 
shown here with the Zheng et al. (1997) UV behaviour, 'UVHST', and 
extended into the IR as set out in ||3;2J The solid line 'RRl' and dotted 
line 'R,R,2' are the enrpirical AGN templates based on Rovi'an— R.obinson 
(1995), shown with either a drop-off in the UV ('UVdrop', dashed) or a rise 
in the UV ('UVrise', solid). The 912A Lyman limit has been indicated, as 
is the slope of a power-law continuum \vith ck_^ — — 1.5. Discussion of these 
AGN is covered in *|^^j2J with UV treatment discussed in J^^ 



The Sab, Sbc, Scd and Sdm templates were extended 
in a similar manner, using the sub-912A part of the spec- 
trum (suitably scaled to give a factor 10 rise in flux across 
the 912A discontinuity) from the S cd template of CWW, ex- 
tended by iBolzonella et al.1 11200(1) . They were all given the 
same UV behaviour because the difference in UV spectrum 
between these galaxy types is expected to be less than the 
uncertainty in their actual UV behaviour. Indeed at redshifts 
where this region enters the optical filters, the dominant ef- 
fect is due to IGM absorption. 

The starburst templa te was ext ended following t he results 
of Bruzual & Char log ||l993) and lFioc &: Rocca-V olmcrangg 
lll997D which show that within lOMyr the UV light drops 
sharply due to the evolution of massive stars off the Main 
Sequence. Below 912A the starburst was assumed to be opti- 
cally thick, an assumption that has been verified for nearby 
starburst galaxies by Loithercr ct al. ( 1995) . 

It is noted here that as well as the UV behaviours 
used above for the 6 galaxy templates, other forms were 
tested, such as sharp cut-offs at 912A - essentially as- 
suming that the galaxies are optically thick to ionizing 
radiation below the Lyman limit - or simply taking the 
flux at lOOOA and setting this value for sub-lOOOA wave- 
lengths. Such approaches were not found to be as successful. 



this is perhaps a more straight-forward task, since the emit- 
ted electromagnetic energy can be assumed to be due to 
stars and dust heated by those stars. There exist a num- 
ber of stellar synthesis codes that, when coupled with spec- 
tral evolution, can self-consistently reproduce the emission 
of galaxies across many orders of magnitude in wavelength. 
The SSPs adopted by the spectral synthesis code and used 
to generate the 6 galaxy templates in this work were not 
detailed enough below lOOOA; hence for the far-UV part of 
these templates the results of other spectral evolution codes 
has been considered. 

B oth the isochrone synth esis of [ Bruzual fc Chariot 
(Tggi) and the PEGASE code of lFioc fc Rocca-Volmerange 
(1997) show that for old stellar populations (older than sev- 
eral Gyr) such as those that characterise ellipticals, there 
is a rise in the far-UV due to low-mass stars in their post- 
AGB evolution. Furthermore, hot post-AGB stars decrease 
the amplitude of the 912A break once their envelopes have 
dissipated. 

In order to extend the elliptical template into the far- 
UV, therefore, the sub-912A part of the spectrum from 
the elhptical template in HYPERZ (BolzonoUa ct al] l2000l) 
was us ed, who had extended it from t he elliptical tem- 
plate of lColeman. Wu fc Weedmanllll980l) . hereafter CWW. 
The flux was scaled in order to give a slight rise across 
the 912A discontinuity. The four template set of CWW 
(elliptical, Sbc, Scd, Im) h as been used in many ph oto- 
metric redshift studi es fE.g lGwvnlll999l fBenitedll999l and 
iBrodwin et al.|[l999|) and has been found to be a robust 
and reasonably complete t emplate set. The extension into 
the far-UV carried out by 'Bolzoncll a et al] i2000l) was by 
means of [Br uzual fc Chario t (1993) spectra with parame- 
ters (SFR and age) selected to match the o bserved spectra 
at zer o redshift. They infact use the IMF of lMiller fc Scalo| 
il979h . however this choice has a negligible impact on the 
final results, as they discuss in their §4.6. 



Determining the exact shape of the UV spectra of AGN 
is problematic for the same reasons as for galaxies. Addition- 
ally, the observed broad continuum feature in the optical- 
UV, the 'big blue bump', is confused by many broad and 
blended lines, which are thought to be due to fast moving 
ionised clouds near the centre of the AGN. Contamination 
from the host galaxy also has an effect for low luminosity 
AGN, particularly in the case of a nuclear starburst. Hence 
several far-UV trends have been investigated to see how they 
effect the accuracy of quasar photometric redshifts: 

(i) although observation s of the cont i nuum blueward of 
the Lyman edge are rare, IZheng et al.l lll997l) constructed 
a composite spectrum from 284 Hubble Space Telescope 
(HST) Faint Object Spectrograph (EOS) spectra of 101 
quasars at 2 « 1. Around 90 per cent of the sample were at 
redshifts < 1.5 so the region blueward of Lya could be stud- 
ied without large effects from the Lya forest. The shortest 
wavelength data, 350-600A, were drawn from higher red- 
shift quasars for which significant corrections for the Lya 
forest and continuum absorption were made. There appears 
to be a break in the continuum slope at around 1050A with 
ax « —0.2 in the far UV. This composite spectra was used 
as one possible UV extension to the sub-912A SEDs of the 
'SDSS', 'RRl' and 'RR2' AGN templates and can be seen 
in Eig.Elas the line labelled 'UVHST'. 

(ii) the 'UVHST' is approximately a flat continuum. In 
order to explore the two alternatives, the sub-912A SEDs of 
the 'SDSS', 'RRl' and 'RR2' AGN templates were also mod- 
elled as: a 'UVrise' with cia = — 1.5; a 'UVdrop' with a\—+3. 
These can also be seen in Eig.|3 Recall that IGM absorption 
is applied to these SEDs depending on redshift. 



4 SPECTROSCOPIC COMPARISONS 

An important stage in the development of any photometric 
redshift code is to run it on a catalogue of objects for which 
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the redshifts are already known from spectrocopic obser- 
vations. A version of the ImpZ code has already used to 
study the evolution of the UV radiation density, the dust 
opacity, and hence the SFH for galaxies in the HDF North 
and South (RR03). As part of this, a comparison was made 
between the photometric output and the spectroscopic red- 
shifts of 152 HDF North galaxies. For galaxies with at least 
4 photometric bands (U,B,V,I) it was found that the spec- 
troscopic redshifts were successfully matched to an accuracy 
of around 10 per cent in (1+z). Around 2.5 per cent (defined 
where there was one or more secondary minima with x^ less 
than 1.0 above the global minimum and (H-z) differing by 
more than 20 per cent) of galaxies were found to have a 
'significant' redshift alias . 

The version of ImpZ in RR03 has also been used to gen- 
erate photometric redshifts for the final band-merged ELAIS 
catalogue (6.7, 15, 90 and 175/im, and U , g' , r' , i', Z, J, H, 
K and 20cm; ?), fitting both galaxies and quasars by using 
their U , g' ,r' ,i' , Z, J, H , and K data. Again the photomet- 
ric redshifts were accurate to around 10 per cent in (H~z), 
with a greater dispersion in Az for AGN fits. 

In order to extend the applicability of ImpZ, this inves- 
tigation explores Impz's application to two further sets of 
photometric data with spectroscopic redshifts. 



the galaxy templates) based on class flags. The INT WAS 
sources are classed in each band as either —1: point-like, 0: 
noise, 1: non-stellar, —9: saturated, —2: could be stellar, —3: 
might just be stellar. The best bands to deflne stellarity are 
g' , r' and i' so if the flag is —1 in any of these bands then 
the source is defined as stellar. The choice of this definition 
of stellar was reached after extensive tests with the spec- 
troscopic INT WAS ELAIS catalogue which found that this 
identified the most rmmber of actual quasars whilst keeping 
galaxy contamination to a minimum. 

This stellar definition splits the sample into 24 stellar 
and 139 non-stellar sources. 

Galactic extinction in ELAIS Nl and N2 is low (they 
were after all selected for their low 100^.m intensity). 
Both have extinctions of E(B — V) = 0.007, which using 

I . I I r 

the extinction-wavelength relation in lCardelli et alj 1119891) 
equates to Av ~ 0.022 for Ri,=3.1. We can also use the same 
relation to find the extinction in each of the WAS bandpasses 
and correct for it (see Table Ol. The effect of varying R„ on 
the shape of the extinction curve is most apparent at the 
shorter wavelengths, and along with uncertainties in the ac- 
tual form of the average R„-dependent extinction law, the 
accuracy of these corrections is around 0.002 in magnitude 
(less than errors in the actual photometry). 



4.1 ELAIS spectroscopic redshifts 

The spectroscopic data in the ELAIS Nl and N2 fields 
(Perez-Fournon et al. 2003; in prep.) is comprised of two 
samples: 

a) objects observed with the William Herschel Tele- 
scope and the fibre-fed WYFFOS/Autofib2 spectrograph at 
the Observatorio del Roque de Los Muchachos, La Palma, 
as part of an International Time Project approved by the 
Comite Cientifico Internacional of the Observatories of the 
Instituto de Astrofisica de Canarias for the year 2000 (PI I. 
Perez-Fournon), and 

b) objects from the SDSS First Data Release 
JAbazaiian et alll2003t) . 

These data comprise 172 extragalactic sources, with 
spectroscopic redshifts ranging from a 0.0264 to 2.9426, and 
detections in to 5 of the INT WAS bands. The mean limit- 
ing magnitudes (Vega, 5a in 600 seconds) in each filter across 
these two areas are 23.40(U), 24.94(g'), 24.04(r'), 23.18(i'), 
and 21.90(Z). The spectroscopic redshift information also 
includes a 'QSO'/'GALAXY' flag (26/146 objects). These 
data are reduced to 163 by removing objects flagged as sat- 
urated via the INT WAS class flags (1 source), objects that 
have no WAS detections (4 sources), objects with only one 
detection (3 sources), and objects flagged as being contam- 
inated by a bright star (1 source), leaving 157 objects with 
detections in 5 bands, and 6 with detections in 4 bands, 
sub-divided into 25 'QSO' and 138 'GALAXY' sources. The 
'QSO'/'GALAXY' information for INT WAS data enables 
us to decide on the best combination of the INT WAS class 
flags for deflning possible AGN's that should be included in 
AGN template-fitting, and additionally gives a direct mea- 
sure of how well the galaxy and AGN templates manage to 
separate out these two populations. Hence a choice can then 
be made of which AGN templates to use, and what treat- 
ment in the UV is most successful. Sources are defined as 
stellar (and so AGN templates are considered in addition to 



4.2 CDFN spectroscopic redshifts 

The iBarger et alJ (|2002) catalogue in the CDFN was also 
used. This comprises an X-ray selected catalo gue (from the 
IMs Chandra observation of the HDF North; .Bra ndt et alJ 
l200lfl of 169 objects with spectroscopic redshifts and broad 
band photometry {B, V, R, I, z' from Subaru/Suprime Cam 
and a notched HK' filter with a central wavelength of 1.8/im 
from the University of Hawaii 2.2m telescope/QUIRC - see 
Fig. for filter transmissions combined with the CCD re- 
sponses). The resolution of the X-ray and optical/near-IR 
observations are similar (around 1") with nearly all true 
counterparts expected to lie within the 2" {5a) radii of the 
X-ray source for sources within 6.5' of the approximate X- 
ray image centre and 3.6" radii (4(t) for sources beyond this 
radius. The la limiting magnitudes are (AB magnitudes): 
B(29.0), F(28.5), i?(29.2), 1(27.6), z'(27.0) and HK'{23.3, 
Vega). As no errors are provided with these measurements, 
photometric errors of 0.05 in the Subaru bands and 0.15 
in HK' have been assumed. The catalogue is reduced to 
161 sources after removal of those flagged as saturated (7 
sources) or contaminated (1 source). In addition, HK' de- 
tections were only used if they were brighter than magni- 
tude 20.0 for data quality reasons. This means there are 105 
sources with detections in 6 bands, 52 with 5 bands and 4 
with four bands. Dropout treatment was applied to the two 
shortest wavebands - B and V. 

Since this sample is X-ray selected, we can expect a 
large proportion to be AGN, or AGN-dominated. Hence al- 
though this sample does not have an exactly analogous set 
of filters to the INT WAS, it is another excellent test-bed 
for the AGN template-fitting. It also has a larger sample 
of high-redshift galaxies which test the ImpZ application 
to galaxies across a broader redshift range. In place of the 
stellar definition in the application to the INT WAS sur- 
vey, AGN template-fitting was carried out in addition to the 
usual galaxy template-fitting for objects that are flagged as 



Table 3. ELAIS Nl and N2 Galactic extinction corrections for tlie INT 
WAS filters and CDFN Galactic extinction corrections for the Bargcr ct al. 
(2002) filter set, derived using Cardelli et al. (1989) with R„=3.1. 



INT 


WAS 


filters 










filter 


U 


g' 


r' 


i' 


Z 




Aa 


0.034 


0.025 


0.019 


0.013 


0.010 




Barger 


et al. 


filters 










filter 


B 


V 


R 


/ 


z' 


HK' 


Aa 


0.049 


0.037 


0.031 


0.022 


0.018 


0.005 



being optically compact [C ~ 23 sources], having broad- line 
features [B - 6 sources], or having both [BC - 23 sources]. 
Hence 52 of 161 sou rces have AGN template-fitting applied. 
iBareer et all ll2002ll used HYPERZ on this sample but only 
managed to get about 1/3 of their broad-line sources' pho- 
tometric redshifts within 25 per cent of the spectroscopic 
values. 

Galactic extinction in GDFN is again low, with ex- 
tinction E(B — 1/)=0.012, which using the extinction- 
wavelength relation in lCardelli et al.l il989fl equates to Av ~ 
0.037 for Ri,=3.1. We can also use the same relation to find 
the extinction in each of the bandpasses and correct for it 
(see Table EJ. 

4.3 Comparison to HYPERZ 

In order to see how the popular template-fitting photomet- 
ric code HYPERZ compares, it has been run on the same 
two catalogues, using as similar parameters as possible. The 
set of 6 galaxy templates was used, with RRlUVrise and 
RR2UVrise used for AGN-fitting. One difference between 
the HYPERZ and ImpZ codes is that the A„ prior that in- 
creases the Xred for increasing ]A„] is not implementable in 
HYPERZ. Similarly the prior that makes AGN fits more 
preferable than galaxy fits for stellar sources is not avail- 
able in HYPERZ (these priors were set out in ^TH . HY- 
PERZ was run with the same redshift range, but a redshift 
step in (1-fz), rather than in logio(l-l-z). The same absolute 
magnitude limits as in ImpZ were used (but the equivalent 
limits in g' band for INT WAS) for the galaxy and AGN 
fits, with an A„ range of 0.0 to 1.0 for galaxy fits and A„=0 
for AGN and elliptical fits. Two reddening laws were tried - 
[Calze tti ct al. (20o3) and|£eaton (1979) (fit by Fitzuatrick 
Il989l for Milky Way), w ith galaxies getting the best results 
with the ISeatonI il979l) law. Allowing A„ to vary for AGN 
fits gave very similar results to setting Au=0, with the best 
results shown in Fig. 2] and Fig. |S| The results of this com- 
parison to HYPERZ are in JO] 



5 SPECTROSCOPIC REDSHIFT RESULTS 

We measure the reliability and accuracy of the photometric 
redshifts via the fractional error Az/(1 + z) for each source, 
examining the mean error Az/(l + 2), the rms scatter Uz and 
the rate of 'catastrophic' outliers r), defined as the fraction 
of the full sample that has ]Az/(l + z)\> 0.2. 
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l^z/{l-\-z), Az / (1 + z) and u^ are calculated as follows: 
Az/{l + z)= hp'-"^ - ^-p-- \ (3) 



and 

Az/il + z) = J2[^ 
and 
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with N being the number of sources with both spectroscopic 
redshifts and photometric redshifts. 

For both spectroscopic samples, the outlier-clipped 
rms, ared, calculated from sources with |Az/(l -1- 2;)] < 0.2 
was around 0.07 and the outlier fraction was 4.9 per cent 
for the INT WAS ELAIS and 12.4 per cent for the CDFN 
sample. 

5.1 Results of ELAIS spectroscopic study 

The code parameters and templates that gave the best re- 
sults are described in t|5.5l Fig.|l]shows the results of running 
ImpZ on the INT WAS ELAIS spectroscopic redshift sample 
with this 'optimum' set of parameters and templates, with 
Fig. 13 showing in close-up the behaviour at low redshift. As 
weU as plotting logio(H-Zspec) versus logio(l+Zphot) for the 
whole sample, the results for just those with a spectroscopic 
'QSO' flag are shown, as is the 'GALAXY' sub-sample. This 
better illustrates the varying success of the code on these dif- 
ferent groups of sources, and the template-types that were 
best-fitting. In particular, it can be seen that 20 of the 25 
'QSO' objects were best-fitted by AGN templates. Four of 
the five that weren't were defined as non-stellar following 
the definition set out in i|2.2l (so AGN-fitting was not applied 
to them) and all five are at low redshift {z < 0.3). 

ImpZ found solutions for all but one of the 162 sources. 
The failed source is discussed in t|5.l.5| 

For the sample as a whole, the total rms scatter, atot, 
was 0.12, with Az/{l + z) = —0.02. The outlier-clipped rms, 
(Jred was 0.068 and the rate of 'catastrophic' outliers, r] was 
4.9 per cent. 

The code was successful in identifying and fitting AGN 
templates to nearly all of the 'QSO' objects, and was quite 
successful at returning an accurate Zphot- The total rms scat- 
ter, atot, was 0.29 for this sub-sample, with Az/{1 -j- z) = 
—0.11. The outlier-clipped rms, ared was 0.11 and the rate 
of 'catastrophic' outliers, rj, was 32 per cent, higher than for 
the spectroscopic sample as a whole. It can be seen from 
Fig. 2] however that infact the agreement with Zspec is rea- 
sonable (all but one are within the Sarsd limits), especially 
considering the infancy of applying photometric redshifts to 
quasars. Nine of the 25 'QSO' sources were fit within 25 per 
cent of their Zspec values. 

The code was even better on the 'GALAXY' sources, 
with atot=0-061, Az/(1 + z) = -0.029 and ared^O.061. 
There were no 'catastrophic' outliers. Although these 
'GALAXY' sources were all relatively low redshift (z < 0.2), 
the success is encouraging. 

Statistical results for both the 'best-case' setup and also 
for other ImpZ setups can be found in Table 2] 
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Figure 4. Spectroscopic validation results for the INT WAS ELAIS and CDFN spectroscopic samples. The rms values shown arc the outlier-clipped rms 
values, fTj,g^. Top two rows shows INT WAS results (first row is ImpZ results, second row is HYPERZ results), set out as follows: first plot is all the data 
(crosses) with those flagged as 'QSO' outlined with a square. The dotted lines are the Zcr^,^^ limits; second plot is only 'QSO' objects (crosses) with those 
fit as AGN outlined with a square. Dotted lines are again SCTj^^d- fouth plot are 'GALAXY' objects (crosses) with those fit as AGN outlined by a square. 
Dotted lines arc 3cr,,ed limits. The third and fourth rows shows results from the CDFN sample with a similar format as to the first row: first plot is all the 
data (crosses) with those flagged as 'B', 'C' or 'BC outlined with a blue square. The dotted lines arc the 3(7^-g^ limits; third plot is only 'B', 'C or 'BC' 
objects (crosses) with those flt as AGN outlined with a square. Dotted lines are again Scr^g,^; fourth plot shows remaining objects (crosses) with dotted lines 
the Z(Jj,^fi limits (third row is ImpZ results, fourth row is HYPERZ results). 
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5.1.1 Stellar definition 

It is important to restrict which sources the AGN templates 
are fit to. If the AGN templates are fit to all sources, instead 
of only fitting to those defined as stellar, then the overall 
accuracy of the ImpZ code on the INT WAS spectroscopic 
sample drops to a"tot=0.16 and r) rises to 6.8 per cent. This 
is because the inclusion of AGN templates introduces more 
degeneracies into the color-redshift space and we therefore 
wish to use further information, in this case the class fiags, 
to break some of these degeneracies. 

Similarly, although the stellar sources are more likely 
to be quasars, if we restrict the fitting of galaxy templates to 
non-stellar objects only, then the accuracy of ImpZ results 
for stellar objects is reduced because not all of them will 
be quasars. If galaxy templates are only fit to non-stellar 
sources then crtot=0. 16, A2/(l ~\- z) — 0.01, crred=0.07 and 
?7=7.5 per cent. Instead, galaxy templates are fit to stellar 
objects, but AGN template-fits are made preferable through 
the prior set out in i l2.2l The alteration of this prior is dis- 
cussed in t|5.1.4l 

Bearing in mind that the code is to be applied to the 
INT WAS ELAIS survey as a whole where there will be pro- 
portionally less quasars and more galaxies, we wish to make 
the criteria for a stellar definition and ensuing AGN-fitting 
as tight as possible in order to minimize galaxy contamina- 
tion (and improve efficiency). The best combination of class 
fiags for identifying quasars without undue contamination 
from galaxies was found to be stellar if the class fiag is —1 
in g' , r' , or i' band. This defines 21/25 'QSO' sources and 
3/138 'GALAXY' sources as stellar. 



5.1.2 AGN fitting 

The choice of which AGN templates to include in the 
template-set, and to a lesser extent how to treat them in 
the UV, had a large effect on the success of the code in 
obtaining accurate redshifts for stellar objects, and also in 
whether the 'QSO' sources were correctly identified as such 
and best-fitted by an AGN rather than a galaxy template. 
At high enough redshifts for the different UV treatments to 
enter the INT WAS bandpasses {z > 2.5) there are relatively 
few sources available to test the different behaviours - the 
CDFN sample is more useful (see tl5.2.H . For the highest 
redshift INT WAS source ('QSO' a.t z = 2.94) it was found 
however that the AGN templates were not the best-fitting 
when the 'UVHST' behaviour was used, with instead a low 
redshift Sbc template being fit. 

The 'UVdrop' behaviour was found to be reasonable, 
with the high-redshift source being fit as an AGN at z = 3.4 
whilst the 'UVrise' gave the redshift as z — 3.1. The 'UVrise' 
behaviour was therefore chosen to be the more successful, al- 
though the main reason for its choice was determined by the 
CDFN results (' «)5.21 which had more high-redshift sources. 

The 'SDSS' AGN template was found to best-fit a rea- 
sonable number of stellar sources but unfortunately the re- 
sulting photometric redshifts were inaccurate, with a larger 
number of catastrophic outliers. For example, using the 
SDSS template with the 'UVrise' behaviour best-fit 19/25 
'QSO' sources, but with an increased a"tot=0. 30 and an out- 
lier fraction of more than half. 

It was found that the best combination of AGN tem- 



plates was to use the RRl and RR2 templates with the 
UVrise behaviour. Using just one or the other tended to in- 
crease the rms scatter of the 'QSO' objects from atot~0.2S 
to roughly 0.40 and to increase the number of outliers. 

It is noted that if no AGN templates are used — only 
the 6 galaxy templates — then the outlier fraction for the 
'QSO' sources increases to ri=80 per cent with crtot=0. 48. 
Indeed, no source with Zspec >0.6 is fit within Satot of Zspec. 



5.1.3 High redshift sources 

Since the main power of photometric redshifts is in their 
application to higher redshifts, it is informative to examine 
the higher redshift sources as a separate sample. Taking our 
redshift cut as Zspec > 0.3 (sources below this redshift are 
plotted in close-up in Fig.|^ means that there are 22 sources 
in the high redshift sample, comprising 2 'GALAXY' and 
20 'QSO' sources. Hence the statistics for this sample are 
similar to that for the QSO sources alone, with atot=0.34, 
Az/{1 + z) = 0.12 and (Jred^O.ll with an outlier rate of 
36 per cent. As would be expected, then, the accuracy for 
higher redshift sources is less than for lower redshift sources, 
but is still good. A better grasp of high redshift performance 
is provided by the CDFN sample, which has more objects 



at Zspec > 0.3 - see >I5.2.2I 



5.1.4 Relaxation of parameters 

5.1.4.1 At, If Av fitting is completely turned off, so 
that only At,=0 solutions are allowed then the majority 
of the increase in rms comes from increased scatter in the 
'GALAXY' sub-sample. Hence the inclusion of Ai, as a pa- 
rameter improves the accuracy of the redshifts, whilst pre- 
sumably giving some information about the extinction of 
each source. Several different At, limits were tried for the 
galaxy templates. For example, since the ELAIS fields were 
originally IR-selected we might expect some sources to have 
high extinction so a range of A„ =0-3 was used, with /3 = 3 
instead of 2 in the A,, prior that minimises Xred + P-^v- This 
did not alter the results.. However, extending At, freedom to 
AGN templates increases the the outlier-rate for 'QSO's to 
V = 40%. 



5.1.4.2 AGN prior If the prior that a galaxy-fit to a 
stellar source must have Xred-""'""" < iXred-^^'^ ~ 4) is 
removed then only 17 of 25 'QSO' sources are fit as AGN and 
the highest redshift source ('QSO' at z — 2.94) is instead fit 
as a low-redshift Sbc galaxy, with an overall increase in the 
rms and outlier fraction. If the prior is reduced in strength, 
to X^d-""'""" < ixled-'^^'^ - 2) then although the overaU 
statistics are almost as good, the highest redshift source is 
again fit as a low-redshift Sbc galaxy. 



5.1.4.3 Absolute Magnitude limits The choice of 
the absolute B-band magnitude, M_g, limits of the photo- 
metric solutions is important for removing physically un- 
likely low-redshift low-luminosity and high-redshift high- 
luminosity solutions which greatly increase the number of 
catastrophic outliers. For example, allowing brighter AGN 
limits of -29<Ms<-17 instead of -27<Ms<-17 resulted 
in a number of high photometric redshift (2 > 5) fits to 
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low redshift sources. Similarly, increasing the galaxy lim- 
its to -25.5<Ms<-13 replaced a number of correct faint 
low-redshift fits with incorrect brighter high-redshift fits. In- 
correct low-redshift solutions for high-redshift souces were 
found to increase if the opposite action was taken and 
brighter galaxy and AGN fits were not allowed. The use of an 
upper Mb envelope that increased with redshift was found 
to be more successful than using a fixed value for galaxies, 
as might be expected from the known strong evolution in 
galaxy luminosities. The selected redshift dependence used 
here was [-22.5 - 2logw{l + z)] <Ms < -13.5, though 
[-22.5 - 2/3] <Mb < -13.5 worked equally well. 
The final absolute magnitude limits were chosen in order to 
give the best agreement with the spectroscopic data whilst 
taking into consideration that the INT WAS ELAIS survey 
as a whole will have a greater diversity of sources. 

5.1.5 Outliers 

ImpZ fails to find a solution for one INT WAS ELAIS source, 
a 'GALAXY' at redshift 0.1454 with a non-stellar fiag. This 
source is not detected in g' , but is relatively bright in other 
bands (19.42, 18.05, 17.18, 16.75 in U, r' , i' and Z). It was 
flagged as having multiple counterparts so this might imply 
that the photometry is incorrect. HYPERZ manages to find 
a solution at Zhyp~0.19 as an Sab with A„=0.4. 

The remaining most conspicious outliers are a group 
of three sources with high spectroscopic redshifts and 
lower photometric redshift solutions, falling below the 3ared 
boundaries. It is noted that two of these sources are at red- 
shifts where the 4000A break has left the longest INT WAS 
waveband (at z ~ 1.4) but the 912A Lyman limit has still 
not entered the shortest waveband {z ~ 2.5). Sources in 
this redshift range are harder to successfully derive redshifts 
for since the primary SED features on which photometric 
redshift-fitting relies are not available. Table Q details the 
redshifts where common spectral features enter and leave 
the bandpasses used in the INT WAS and CDFN catalogues. 
For five other sources that fall into this redshift range, ImpZ 
is still successful in deriving the correct redshift, implying a 
success-rate in this redshift region of 71 per cent. 

The other outlier (a 'QSO') lies above the 3ared bound- 
aries, with a low Zspec and a high Zphot solution as an AGN. A 
solution near the correct redshift is found if the SDSS AGN 
template is used, however including this template degrades 
the performance of the code for the rest of the sample. 

5.2 Results of CDFN spectroscopic study 

Fig.|l]shows the results of running ImpZ on the CDFN spec- 
troscopic redshift sample with the same 'optimum' set of 
parameters and templates as plotted for the INT WAS sam- 
ple. As expected, since there is less available information 
on which to base the stellar and non-stellar definition and 
no actual 'QSO' or 'GALAXY' flag, the resulting analysis 
can be less clearly separated into results for galaxies and 
quasars. 

ImpZ found solutions for all 161 sources. The total rms 
scatter, atot, was 0.17, with Az/{l+z) = —0.01. The outlier- 
clipped rms, Ored was 0.07 and the rate of 'catastrophic' 
outliers, 77 was 12.4 per cent. From Fig.^Jit can be seen that 
the code really only went badly wrong for 3 sources. 




Spectroscopic Redshift 



Figure 5. Zoomed view of spectroscopic validation results for the INT 
WAS ELAIS spectroscopic sample. First plot sho\vs all ImpZ data for 
2spec < 0.3 (crosses) w^ith those flagged as 'QSO' outlined ^vith a square; 
second plot is the same plot for the HYPERZ results. For both plots the 
scatter appears large since the photometric redshift method is accurate to 
perhaps 0.05 in (1+z) — of order of (1+z) for low rcdshfits. 



The code was successful in identifying and fitting AGN 
templates to 27/52 of the broad-line or compact objects, and 
overall was successful at returning an accurate photometric 
redshift. The total rms scatter, otot, was 0.19 for this sub- 
sample, with Az/(1 + z) = —0.01. The outlier-clipped rms, 
(Jred was 0.07 and the rate of 'catastrophic' outliers, 77, was 
17.3 per cent. 

The code is also successful for the remaining sources 
(not B and/or C sources), with atot =0.16, Az/(1 ~\- z) = 
—0.01 and (7red=0.07. The rate of 'catastrophic' outliers, 77, 
was 10.1 per cent. 

Statistical results for both the 'best-case' setup and also 
for other ImpZ setups can be found in Tabled 



5.2.1 AGN fitting 

Nineteen of 29 (66 per cent) broad-line (B or BC) source 
redshifts were within 2 5 per cent of t heir spectroscopic val- 
ues. As a comparison, iBarger et alJ ll2002i'l got 1/3 of the 
sample within this tolerance. 

If we use just BC (broad-line and compact) or C (com- 
pact) as the stellar definition (in plac e of B, C, BC) then 
the broad- line source numbered 174 in lBarger et al.l 112002^ 
(with the highest spectroscopic redshift, z = 5.19) is incor- 
rectly fit as a low-redshift galaxy. The reason why the B 
flag should perhaps not be used in the stellar definition is 
that this is a spectroscopically rather than photometrically 
derived property, and so would not be available to a photo- 
metric survey in general. 

As with the INT WAS ELAIS spectroscopic sample, 
the most successful AGN templates were found to be the 
RRl and RR2 templates and due to the greater number 
of high-redshift sources in the CDFN sample, the best UV 
treatment could investigated more clearly. It was found that 
using the UVdrop or UVHST behaviours was less successful 
than UVrise. The UVHST behaviour increased the popula- 
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Figure 6. Spectroscopic validation results for the INT WAS ELAIS and CDFN spectroscopic samples. Top two rows show INT WAS results (first row 
shows ImpZ results, second row shows HYPERZ results), set out as follows: first plot is the distribution of residuals betiveen Zgpec and Zp^^^, normalised to 
(l+Zspec)- Solid histogram is for the wrhole sample, dashed histogram is that of QSO objects only; second plot shows the reduced chi distribution of the 
fits; third plot is the type distribution of the fits where 1 = E, 2 = Sab, 3=Sbc, 4=Scd, 5 = Sdm, 6=Sb, 7=RRlUVrise and 8=RR2UVrise; fourth plot is the A„ 
distribution. The third and fourth rows show results from the CDFN sample with the same format as the first tw^o rows except that the dashed histogram is 
now that of broad-line and/or compact sources only (third row shows ImpZ results, fourth row shows HYPERZ results). 



tion of low-redshift sources incorrectly fit as high-redshift, 
although the overall success was good, with the main result 
being a slight increase in the outlier fraction for the stellar 
sources. The UVdrop treatment had a similar effect, making 
the UVrise behaviour the most effective. 

It is noted that if no AGN templates are included in 
the template-set then the outlier-rate for stellar sources 
increases to 58 per cent, with a large population of high- 
redshift stellar sources incorrectly placed at lower redshifts. 
If AGN templates are fit along with galaxy templates to non- 
stellar sources then the outlier-rate for non-stellar sources 
increases to 28 per cent and atot rises to 0.88. If, conversely, 
galaxy templates are only fit to non-stellar sources then the 
rms scatter of the stellar sources becomes much larger, with 
atot=lA2. 



5.2.2 High redshift sources 

Again, we examine the higher redshift sources as a separate 
sample. Taking our redshift cut as Zspec > 0.3 (sources below 
this redshift are plotted in close-up in Fig. (SJ means that 
there are 143 sources in the high redshift sample, comprising 
92 non-stellar and 51 stellar sources. The statistics for this 
sample are atot=0.18, Az/{1 + z) — 0.01 and (Tred=0-07 
with an outlier rate of 13 per cent. The accuracy for higher 
redshift sources is therefore still good. The performance is 
better than for the high-redshift sample in the INT WAS 
ELAIS catalogue since in that case 90 per cent of the sources 
were 'QSO', for which photometric redshift-fitting is less 
accurate. 



5.2.3 Relaxation of parameters 

5.2.3.1 Ai, If A^ fitting is completely turned off, so that 
only A„=0 solutions are allowed then rj increases to 17.4 per 
cent. Again, the inclusion of A„ as a parameter improves the 
accuracy of the redshifts, though analysis of a sample with 
known Ai, would be required to quantify how the resulting 
A„ of the solution compares to actual A„. 



5.2.3.2 AGN prior If the prior that a successful galaxy- 
fit to a stellar source must have xLd-^"'"''^ < (Xred-^'^^-4) 
is removed then only 13 of the 52 stellar sources are fit 
by AGN and many of those that are then fit as galaxies 
are placed at much lower Zphot- Using a reduced prior of 
^2^^ saiaxB ^ (^2^^ AGiv _2) ^auses the highest Zspec source 
to again be fit as a low redshift galaxy, and increasing the 



prior to Xr 



galaxy 



< (X,e 



25) creates a number of 



false high-Zphot AGN solutions. 



5.2.3.3 Absolute Magnitude limits As with the INT 
WAS sample, allowing brighter AGN limits of -29<Mb<- 
17.5 instead of -27<M_b <-17.5 resulted in a number of high 
photometric redshift (2 > 5) fits to low redshift sources. Sim- 
ilarly, increasing the luminous galaxy limits to -25.5<M_b <- 
13.5 replaced a number of correct faint low-redshift fits with 
incorrect luminous high-redshift fits. Incorrect low-redshift 
solutions for high-redshift souces were found to increase if 
the opposite action was taken and brighter galaxy and AGN 
fits were not allowed. Again, allowing an increase in the up- 
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Table 4. Statistical results for application of the ImpZ code to the INT WAS ELAIS catalogue, for different ImpZ setups. 



INT WAS ELAIS Statistics for different ImpZ setups 


Az/{l + z) 


All sources 


v(.%) 


'QSO' sources 'GALAXY' sources 
Az/(l + z) atat CT^ed V{%) Az/(l + z) 


o-tot 


O-red 


V (%) 


-0.02 


0.12 


0.068 


4.9 


'Best-case' ImpZ setup as in §5.1 
-0.11 0.29 0.11 32.0 


-0.029 


0.061 


0.061 


0.0 


-0.01 


0.16 


0.07 


6.8 


AGN fit to all sources as in §5.1.1 
-0.05 0.36 0.11 40.0 


-0.002 


0.07 


0.06 


0.7 


0.01 


0.16 


0.07 


7.5 


Galaxies only fit to non-stellar sources as in §5.1.1 
-0.08 0.32 0.11 36.0 


0.03 


0.11 


0.06 


2.2 


-0.04 


0.13 


0.07 


SDSS t 
8.0 


cmplate with UVrise as the only AGN template, as in 
-0.21 0.30 0.11 52.0 


1 §5.1.2 
-0.004 


0.06 


0.06 


0.0 


-0.01 


0.16 


0.06 


RRl template with UVrise as tile only AGN template, as in 
6.2 -0.07 0.39 0.10 40.0 


§5.1.2 
-0.004 


0.06 


0.06 


0.0 


-0.03 


0.15 


0.07 


RR2 template with UVrise as the only AGN template, as in 
8.0 -0.17 0.36 0.12 52.0 


§5.1.2 
-0.004 


0.06 


0.06 


0.0 


-0.02 


0.12 


0.07 


4.9 


UV drop behaviour in place of UVrise, as in §5.1.2 
-0.11 0.29 0.12 32.0 


-0.004 


0.06 


0.06 


0.0 


-0.02 


Th 
0.13 


c 4 CWW templates in place of the E, Sab, Sbc, Scd, and Sdm tem 
0.07 5.6 -0.12 0.28 0.12 32.0 


plates, as in 
-0.01 


§6.2 
0.07 


0.07 


0.7 


-0.05 


0.14 


0.07 


Sab, 1 
8.6 


RlUVrise and RR2UVrise as the only templates, as ir 
-0.19 0.30 0.12 48.0 


L §6.2 
-0.03 


0.07 


0.06 


1.5 


-0.07 


0.20 


0.06 


12.3 


No AGN templates, as in §5.1.2 
-0.42 0.48 0.11 80.0 


-0.004 


0.06 


0.06 


0.0 


-0.03 


0.13 


0.07 


6.8 


Ai; free for AGN templates 
-0.17 0.30 0.11 44.0 


-0.004 


0.06 


0.06 


0.0 


-0.02 


0.12 


0.07 


4.9 


A„=0 to 3, with /3 = 3 in A„ prior (see JT^ 
-0.11 0.28 0.11 32.0 


-0.01 


0.06 


0.06 


0.0 


-0.02 


0.13 


0.07 


6.2 


A„ = -0.4 to 1.0 
-0.11 0.29 0.12 36.0 


-0.004 


0.06 


0.06 


0.7 


0.01 


0.13 


0.08 


5.6 


Atj set to zero for all templates, as in §5.1.3 
-0.11 0.29 0.10 32.0 


0.02 


0.08 


0.07 


0.7 


-0.01 


0.15 


0.07 


7.4 


AGN prior removed, as in §5.1.3 
-0.17 0.34 0.12 44.0 


0.02 


0.07 


0.06 


0.7 


-0.01 


0.14 


0.07 


6.2 


AGN prior set to 2 instead of 4, as in §5.1.3 
-0.14 0.32 0.11 36.0 


0.02 


0.07 


0.06 


0.7 


0.003 


0.14 


0.07 


6.2 


AGN prior set to 25 instead of 4, as in §5.1.3 
-0.08 0.32 0.11 36.0 


0.02 


0.07 


0.06 


0.7 


0.01 


0.43 


0.06 


13.0 


Without U band information, as in §6.2 
-0.11 0.43 0.11 60.0 


0.03 


0.43 


0.06 


4.4 


-0.03 


0.13 


0.07 


7.4 


Without Z band information, as in §6.2 
-0.13 0.30 0.10 48.0 


-0.01 


0.06 


0.06 


0.0 


0.01 


0.25 


0.07 


6.2 


Brighter AGN limits of -29<Mb <-17.5, as in §6.3 
0.10 0.64 0.12 40.0 


-0.004 


0.06 


0.06 


0.0 


-0.02 


0.12 


0.07 


Reducing faint AGN limits to -27.0<Mb <-15, as in §5. 
4.9 -O.ll 0.28 O.ll 32.0 


1.3 

-0.01 


0.06 


0.06 


0.0 


-0.02 


0.13 


0.07 


B 

6.2 


righter galaxy limits to -23.5<Ms <-13.5, as in §5.1.3 

-0.11 0.28 0.11 32.0 -0.003 


0.06 


0.06 


1.5 


0.01 


0.33 


0.06 


B 

7.4 


righter galaxy limits to — 25.5<Ms <-13.5, as in §5.1.3 

-0.10 0.29 0.12 36.0 0.03 


0.34 


0.05 


2.2 


-0.02 


0.12 


0.06 


Red 
4.9 


ucing faint galaxy limits to -23<M^ <— 17.5, as in §5 
-0.11 0.28 0.11 32.0 


.1.3 
0.001 


0.06 


0.06 


0.0 



per galaxy Ms envelope with redshift was found to be the 5.2.4 Outliers 

best solution (i.e —22.5 — 2logio{l + z)). 

There are a group of three high spectroscopic redshift 
sources whose photometric solutions are significantly lower 
than their spectroscopic redshifts (there are a total of six 
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sources below the Sared boundary, and one above, but the 
others are close to the boundaries). It is noted that this 
group of three are at redshifts where the 4000 A break has left 
the Z waveband (at z ~ 1.5) but the 912A Lyman limit has 
still not entered the shortest waveband {z ~ 3.3). Sources 
in this redshift are harder to successfully derive redshifts 
for since the primary SED features on which photometric 
redshift-fitting relies are not available. For 20 other sources 
that fall into this redshift range, ImpZ is still successful in 
deriving the correct redshift, implying a success-rate in this 
redshift region of 87 per cent. The slightly better perfor- 
mance than for the INT WAS ELAIS sample (71 per cent) 
for the corresponding 'feature desert' is likely to be due to 
the inclusion of HK' information (which 4000 A enters at 
z « 2.7) and partly due to the poor-number statistics for 
the INT WAS ELAIS sample (only 6 sources in that red- 
shift range). 

5.3 Statistical properties 

The distribution of residuals of the ImpZ solutions in Fig. 
|S| is quite strongly peaked around zero, with a Gaussian- 
like distribution and a slightly non-Gaussian extension to 
larger values. The distributions show that the ImpZ code 
is successful to an accuracy of perhaps 0.1 to 0.2 in (l-l-z) 
for nearly all sources. The distributions of ImpZ solutions 
(dashed line in|SJl for the 'QSO' sub-sample is also Gaussian- 
like, but these distributions are wider and shallower. For 
the INT WAS ELAIS 'QSO's the distribution is peaked at 
around —0.2 rather than zero, implying that the redshift 
solutions tend to be slightly larger than the true values. 
This could be because A„ is set to zero for AGN so a QSO 
reddened by some extinction is fit as unreddened, but at 
a slightly higher redshift. However, this is not seen in the 
larger CDFN sample. 

The Xred distributions of the ImpZ solutions shows that 
the majority of solutions are good, with xi^j^ peaking below 
2. There is no large population of poor-x^ed solutions. From 
the ImpZ type distributions it can be see that Sab solutions 
dominate, with a sizeable fraction of Scd and starbursts and 
very few Sdm solutions. 

The A„ distributions are strongly peaked at zero with 
an even spread of higher values (but no clustering at the 
highest value of A„ = 1 which would imply that the A„ 
range was insufficient). 

5.4 HYPERZ results 

From Fig.0|and Fig.jSjit can be seen that for both the INT 
WAS and CDFN spectroscopic samples, the ImpZ code man- 
ages to do better. The strongest cause of these differences 
are probably the AGN and A^ priors that are used in the 
ImpZ code. In particular, HYPERZ is not as successful for 
the broad-line and/or compact sources in the CDFN sample, 
with a larger population of high-redshift sources incorrectly 
placed at low- redshift (see Fig.^J- 

Overall^ HYPERZ has, for the INT WAS sample, 
atot=0.17, Az/{1 + z) = 0.02, a^ed=0.073 and rj=8.0 per 
cent. For the CDFN sample the statistics were atot~0.26, 
Az/(1 + z) = -0.04, ared=0.Q8 and r?=24.2 per cent. Only 
9/29 broad-line source redshifts were within 25 per cent of 
their spectroscopic values. 



From Fig. |S|it can be seen that the non-Gaussian tails 
of the redshift residuals are more pronounced, particularly 
for the CDFN sample, and the Xred of the solutions are only 
weakly peaked to low values. 

A consideration of the 'feature desert' in redshift space 
where the 4000A break has left the optical wavebands, and 
the 912A Lyman limit has not yet entered shows a marked 
difference from the ImpZ results for the CDFN, indeed it 
is this region where HYPERZ most strongly diverges from 
ImpZ in its outlier-rate: for the CDFN sample, ImpZ had a 
87 per cent success-rate with the 23 sources in this range, 
whereas HYPERZ fails for 10 (57 per cent success). 



5.5 ImpZ summary parameters 

The final best-case setup for application of ImpZ to the INT 
WAS ELAIS as a whole is as follows: 

Templates - 6 galaxy templates (El, Sab, Sbc, Scd, 
Sdm and Starburst) and 2 AGN templates (RRlUVrise and 
RR2UVrise), with IGM and Galactic extinction corrections. 

Fitting - AGN templates fit to stellar sources only and 
galaxy templates fit to all sources, with the prior that a 
successful galaxy fit must have Xred-^"'"''^ < ixled-'^^'^ ~'i) 
if the source is stellar. A„ freedom is allowed for Sab, Sbc, 
Scd, Sdm and Starburst templates, with a prior that makes 
low A„ fits preferable. 

Limits - A„ limits of 0.0 to 1.0 and absolute magnitude 
limits of [~22.5~2logio{l + z)] <M_b<-13.5 for galaxies and 
-27.0<Mfl<-17.5 for AGN. 



6 ERROR ANALYSIS 

Errors in the derived photometric redshifts arise from several 
causes. The first is the inherent error in the measured fiux of 
sources, which is expected to be increasingly important for 
fainter sources. The template-fitting technique adds further 
error due to fitting the continuum of observed galaxy SEDs 
with a set of standard templates which can only sample this 
continuum. This is known as cosmic variance. 

Catastrophic errors are usually due to degeneracies in 
the colour-redshift-extinction space. For example, a late- 
type, heavily extincted galaxy may appear similar to an 
early-type unreddened galaxy, but have a very different red- 
shift. The relative magnitude of these sources of error needs 
to be more properly quantified. 

6.1 Photometry 

In order to remove the effect of cosmic variance from the 
analysis, a synthetic catalogue of fiuxes (L'', g' , r' , i' and 
Z) was generated 5 times, for 7,000 hypothetical sources. 
These sources were created by redshifting the templates 
(equal numbers of each type) and calculating their fiuxes in 
each band. The redshift distribution was chosen to resem- 
ble that of 'main' galaxies in the SDSS First Data Release, 
with a simple Gaussian distribution peaked at —20 with a 
cr of 7 to model the absolute magnitude (r') distribution 
of the sources (additionally, sources were prevented from 
having extreme absolute magnitudes). The fiuxes were ran- 
domly altered in a Gaussian distribution of errors of width 
determined by the la values of the photometric errors in the 
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Table 5. Statistical results for application of the ImpZ code to the CDFN catalogue, for different ImpZ setups. 



CDFN Statistics for different ImpZ setups 


Az/(1 + z) 


All sources 


v{%) 


Broad-line and/or Compact sources 
Az/(1 + Z) atat <T,ed V{%) 


Az/{l + z) 


Other sources 


V (%) 


-0.01 


0.17 


0.07 


12.4 


'Best-case' ImpZ setup as in §5.2 
0.01 0.19 0.07 17.3 


-0.01 


0.16 


0.07 


10.1 


0.17 


0.70 


0.07 


18.6 


AGN fit to all sources as in §5.2.1 
0.01 0.19 0.07 17.3 


0.24 


0.84 


0.07 


17.3 


0.25 


0.80 


0.07 


Galaxies only fit to non-stellar sources as in §5.2 
24.8 0.86 1.42 0.07 56.3 


1.1 

-0.01 


0.16 


0.07 


11.0 


-0.07 


0.23 


0.07 


SDSS template with UVrise as the only AGN template, a! 
22.4 -0.19 0.34 0.08 48.1 


i in §5.2.1 
-0.01 


0.16 


0.07 


10.1 


-0.01 


0.18 


0.07 


RRl template with UVrise as the only AGN template, as 
13.7 0.004 0.21 0.07 21.2 


in §5.2.1 
-0.01 


0.16 


0.07 


10.1 


-0.01 


0.17 


0.07 


RR2 template with UVrise as the only AGN template, as 
12.4 0.01 0.19 0.08 17.3 


in §5.2.1 
-0.01 


0.16 


0.07 


10.1 


0.06 


0.52 


0.07 


14.3 


UV drop behaviour in place of UVrise, as in §5.2. 
0.21 0.88 0.08 23.1 


1 

-0.01 


0.16 


0.07 


10.1 


0.01 


0.36 


The 4 CWW templates replacing E, Sab, Sbc, Scd, and Sdm tem 
0.07 14.3 -0.05 0.20 0.08 21.1 


plates, as in 
0.04 


§6.2 
0.42 


0.06 


11.0 


0.08 


0.57 


0.07 


Sab, RlUVrise and RR2UVrise as the only templates, at 
36.6 0.39 0.96 0.06 59.6 


! in §6.2 
-0.07 


0.21 


0.08 


25.7 


-0.07 


0.31 


0.07 


21.1 


No AGN templates, as in §5.2.1 
-0.18 0.49 0.08 44.2 


-0.01 


0.16 


0.07 


10.1 


0.08 


0.50 


0.06 


23.0 


A„ free for AGN templates, as in §5.2.2 
0.25 0.86 0.05 50.0 


-0.01 


0.16 


0.07 


10.1 


-0.02 


0.18 


0.08 


13.7 


A„=0 to 3, with /3 = 3 in A„ prior {see JT^ 
0.01 0.20 0.08 17.3 


-0.03 


0.16 


0.07 


11.9 


0.03 


0.37 


0.07 


13.7 


A„ = -0.4 to 1.0 
-0.01 0.20 0.07 17.3 


0.05 


0.43 


0.07 


11.9 


0.08 


0.43 


0.08 


17.4 


Ajj set to zero for all templates, as in §5.2.2 
0.22 0.72 0.09 26.9 


0.01 


0.16 


0.08 


12.8 


-0.03 


0.19 


0.07 


13.7 


AGN prior removed, as in §5.2.2 
-0.08 0.23 0.09 19.2 


-0.01 


0.16 


0.07 


11.0 


-0.02 


0.19 


0.07 


14.3 


AGN prior set to 2 instead of 4, as in §5.1.3 
-0.03 0.23 0.08 21.1 


-0.01 


0.16 


0.07 


11.0 


0.06 


0.45 


0.07 


16.8 


AGN prior set to 25 instead of 4, as in §5.1.3 
0.21 0.76 0.08 28.8 


-0.01 


0.16 


0.07 


11.0 


0.13 


0.68 


0.07 


Brighter AGN limits of -29<Mb <-17.5, as in §5.2.2 
14.9 0.43 1.17 0.07 25.0 -0.01 


0.16 


0.07 


10.1 


-0.01 


0.19 


0.07 


Red 

13.7 


Qcing faint AGN limits to -27.0<Ms <-15, as in 
-0.02 0.24 0.08 21.2 


§5.2.2 

-0.01 


0.16 


0.07 


10.1 


0.04 


0.37 


0.07 


Brighter galaxy limits to -23.5<Mb <-13.5, as in §5.2.2 
13.0 0.02 0.17 0.07 15.4 0.04 


0.44 


0.07 


15.4 


0.12 


0.53 


0.07 


Brighter galaxy limits to -25.5<Mb <-13.5, as in §5.2.2 
18.6 0.05 0.20 0.07 19.2 0.16 


0.63 


0.07 


18.3 


0.04 


0.35 


0.07 


Reducing faint galaxy limits to — 23<M^ <— 17.5, as in 
12.4 0.02 0.19 0.08 15.4 


§5.2.2 

0.15 


0.41 


0.07 


11.0 



INT WAS ELAIS catalogue plus 1 per cent of the flux. This 
Monte Carlo simulation generates 35,000 different sources. 
The photometric redshift code is then used to predict, for 
A„=0, the redshifts and spectral types of the sources. Since 
the synthetic catalogue is constructed from the same set of 
templates as those used in the ImpZ code, cosmic variance 
is eliminated, leaving only the effect of photometric error. 

Fig. |7|shows the result of this procedure. The first plot 
(both rows) shows the distribution of residuals between Zi„ 



and Zout, with the first row showing results for low-redshift 
input sources {zin < 1) and the second row showing results 
for high-redshift input sources {zin > 1). The distribution 
of error in redshift for the low-Zi„ sample is very strongly 
peaked at zero and is relatively symmetric about this, but 
the low lying extension to larger errors is not that of a purely 
Gaussian distribution. These large shifts in the redshift oc- 
cur when a secondary minimum in the Xred distribution for 
the source fluxes becomes the primary minimum when the 
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source fluxes are perturbed. Hence for the low-Zi„ sample 
it can be seen that for the large majority of sources, the 
redshift solution is close to the correct value, with only a 
low-level Gaussian spread due to photometric error. 

For the high-Zi„ sample the distribution is also quite 
symmetrical about zero, but is less strongly peaked. This is 
to be expected due to the poorer signal-to-noise photometry 
of these more distant sources. For the high-Zi„ sample, then, 
the dominant effect is the error in the photometry, causing 
the general broadening in the residual distribution. 

The spectral types of the simulated galaxies are com- 
pared with the predicted types in Table|5] This gives the per- 
centage of galaxies for which their spectral type in the syn- 
thetic catalogue was successfully reproduced for each tem- 
plate. 

From this analysis we conclude that without A„ freedom, 
the simulated sources' redshifts are well reproduced within 
0.1 in (zin—Zout) for < z < 1 and roughly 0.4 in (zin-Zout) 
for z > 1. Also, for low redshift sources, photometric errors 
have a reduced effect on the results, whereas for high red- 
shift sources the precision of the photometry has the major 
influence on the redshift accuracy. 

6.2 Templates 

In order to investigate the effect of the choice of template- 
set on the photometric redshifts the 4 UV-extended CWW 
templates (E, Sbc, Scd, Im) are used in place of 5 galaxy 
templates in the ImpZ code (El, Sab, Sbc, Scd, Sdm), re- 
taining the starburst template and the two AGN templates, 
RRlUVrise and RR2UVrise. The code is then applied to 
INT WAS ELAIS and CDFN spectroscopic catalogues. The 
resulting photometric redshifts are compared to the out- 
put when the original templates were used. The results are 
almost as good with atot~0.13, ared~0.07 and 77=5.6 per 
cent for the INT WAS sample and atot=0.36, o-^ed=0.07 
and 77=14.3 per cent for the CDFN. The 6 galaxy templates 
were also replaced with just the Sab template, in order to 
represent an undersampling of SED space. It was found that 
even in this rather extreme example, crtot =0.14, (Jred=0.07 
and 77=8.6 per cent for the INT WAS sample and atot~0-57, 
ared=0.07 and 77=86.6 per cent for the CDFN. 

It is clear that the template-fitting procedure is largely 
unaffected by the specific choice of templates to use, pro- 
vided that the main SED features of normal galaxies (such 
as the Balmer break) are represented. 

Table |S| shows the comparison between the input 
source's template and the template that was fit for it for the 
35,000 sources generated in H6.ll It is immediately clear that 
the majority (75 per cent) of the input template types are 
recovered, and that there is almost no degeneracy between 
early and late-type templates. There is also a clear demarca- 
tion between the galaxy templates and the AGN templates 
- very few AGN templates are fit as galaxy templates (with 
some overlap with the starburst template). Recall that the 
reverse does not occur (galaxy template fit by AGN) since 
sources with an input AGN template are treated as stellar 
by ImpZ. Although there are some degeneracies between the 
templates, this does not effect the photometric redshift ac- 
curacy greatly (as can be seen by the redshift residuals in 
Fig. |7J . Hence the photometric redshift of a galaxy is more 
certain than its exact spectral type, as expected if the bulk 



Table 6. Comparison bet^veen the input model template types and the 
best-fitting templates when ImpZ was run on the Monte— Carol generated 
model catalogues. Values are given as the percentage for each input tem- 
plate that v/as fit by each output template. Ideally only the diagonal region 
(bold) would be populated. 



Out 








Input SED type 










E 


Sab 


Sbc 


Scd 


Sdm 


Sbrst 


RRl 


RR2 


E 


88 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Sab 


8.0 


71.7 


18.6 


0.0 


2.9 


0.0 


0.1 


0.3 


Sbc 


2.0 


8.3 


71.2 


0.0 


0.0 


0.0 


0.0 


0.0 


Scd 


2.0 


8.3 


10.2 


84.8 


4.3 


7.4 


0.0 


0.0 


Sdm 


0.0 


10.0 


0.0 


3.0 


87 


25.0 


0.0 


0.0 


Sbrst 


0.0 


1.7 


0.0 


12.2 


5.8 


67.6 


3.9 


4.7 


RRl 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


64.3 


30.3 


RR2 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


31.7 


64.7 



of redshift identification is due to common features such as 
the Balmer break and in agreement with the results of fitting 
the INT WAS ELAIS catalogue w ith only one gala x y tem - 
plate, and also from the study of iBolzonella et alJ i200Ci) . 



6.3 Filters 

What are the relative importance of different filters? In or- 
der to successfully derive redshifts for different sources, the 
filter set in use needs to have sufficient wavelength coverage 
to encompass the main template features for a broad range 
of redshifts. As seen in H5.1.5l and H5.2.4I most outliers were 
sources whose spectroscopic redshifts meant that the main 
template features (Balmer break and Lyman limit) did not 
fall within filter bandpasses. It is of interest to see how re- 
sults change when a filter is not included in the template- 
fitting. To this end, ImpZ was run on the INT WAS ELAIS 
catalogue without U band information. The statistical re- 
sults of this can be found in Table 2] but the main effect was 
that without U band photometry, all Zphot values were be- 
low 1.5, although 7 sources lie at Zapec > 1.5. This is because 
information on the position of the Lya line has been lost, 
and at these Zspec redshifts no other main template feature 
faUs into the INT WAS bands. 

Additionally, a synthetic catalogue of 5000 sources was 
generated as in >I6.1I but without U band photometry. For 
the Zin < 1 sample, the {zin — Zout) residual distribution 
was strongly peaked at zero, in a similar fashion to the first 
plot in Fig. |7| but with a population of outliers with [zin — 
Zout) ~ —2.6. This suggests that without the use of the 
U band to decide if the Lyman break feature is present or 
not, the Balmer break can be mistaken to be Lya in the 
g' band. Hence the low-Zi„ source is then placed at a high- 
Zout- As a whole, though, the loss of U band information 
does not adversely effect the Zi„ <1 sample. In contrast, the 
Zin > 1 sample is highly dependent on it for a reasonable 
Zout value. The (zin — Zout) residual distribution is almost flat 
- the Zout accuracy is very poor. Thus the U band is central 
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Figure T . Results of Monte Carlo synthetic catalogue investigations with 35,000 sources (Section|^ Top row shows results for 'low' input redshift sources 
(those with z^^i <1): set out as follows: first plot is the distribution of residuals between z^^ and Zo^xt), for template-fitting by ImpZ ^vith A^ freedom turned 
off; second plot is the same, but for solutions with A,; freedom allowed in ImpZ fitting, as described in|5;5J third plot shows the A^ residual from fitting with 
At; freedom turned on. This is simply the A^; of the fit since the input synthetic catalogue was generated without A^, . The second rows displays the same 
information for the 'high' input redshift sources (those with Zj,j >!)■ 



to good photometric redshift accuacy for z > 1.5. This is 
clear from Table Q- without U band none of the template 
features listed is present in the remaining INT WAS filters 
until z > 2.4 when Lya enters the g' band. 



6. 3. 1 Degeneracies in parameter space 

As seen in Hfci.ll and tlfci. 21 there do exist degeneracies in the 
parameter space and this is a cause of outliers in the photo- 
metric error analysis of the simulated catalogues. Common 
degeneracies include confusion between the main template 
features (when either the redshift or the lack of bands means 
most features don't fall into the filter bandpasses) and con- 
fusion between different galaxy templates. Often, however, 
the spectral types can be degenerate without altering z^^ot 
greatly (compare Table |S| with Fig. itTjI. 

It is also of interest to see if A„ introduces its own de- 
generacies. Such degeneracies will be inherent even if A,, is 
not explicitly fit, since real galaxies will have different A„. 
In order to quantify this, the same synthetic catalogue pro- 
cedure as in tl5.1l is carried out, but the ImpZ code is allowed 
to fit for free A„. Since the input catalogue is created with 
Ai,=0 we can see how fitting for free A,, alters the findings. 

Fig. |7] shows the result of this procedure as the second 
and third plots on each row (again the first row is for the 
Zin < 1 sample and the second row is for the Zin > 1 sam- 
ple). The distribution of residuals between Zi„ and Zout are 
very similar to those for the A„=0 case, with both the low- 
zin and high-Zi„ samples being strongly peaked at zero. The 
low-Zi„ sample has a noticeable but low-level tail to higher 
residuals, whilst the high-Zi„ sample can be well represented 
by a broad Gaussian. The high-Zi„ sample has a slight sec- 



ondary bump at {zi„ — Zout) ~ —0.4, which is likely to be 
due to confusion between the Balmer break and 0[III] dou- 
blet, since, as can be seen in Tabled, these features enter 
and leave the i' and Z band with a separation of around 0.4 
in redshift. 

The A„ residual is also shown. This is simply the A„ 
of the fit since the input synthetic catalogue was generated 
without A„. The low-Zi„ sample is very strongly peaked at 
zero, and the high-Zi,i sample is also strongly peaked at zero, 
with slightly more A„ > solutions. This means that the 
increasing photometric error for more distant/faint sources 
causes a few sources to be fit as redder sources than they 
really are. This however does not necessarily cause an in- 
crease in the redshift error, as can be seen from the redshift 
residual plots for the A^-free fitting. It is clear that the ad- 
dition of A„ freedom introduces little degeneracy, and the 
fact that the redshift solutions tend to improve lends weight 
to the argument that the Ai, thus determined reflects the 
physical A„ of the galaxy. 

Running the Monte Carlo simulations with larger pho- 
tometric errors tended to increase the width of the residual 
distributions. 



6.4 Overall errors 

The main sources of error are due to photometry, cos- 
mic variance and degeneracies in colour-redshift-extinction 
space. The relative effect of each depends on the signal- 
to-noise of the source with photometric errors dominating 
those with low signal-to-noise. Splitting into low (2 < 1.5) 
and high redshift {z > 1.5), we can broadly say that the 



ImpZ: a new photometric redshift code 19 



error for low-redshifts is 0.1 or less in (1+z) whereas for the 
high-redshifts the error can be as large as 0.3 or more. 



7 DISCUSSION 

A comparison can be made to the photometric accuracies 
of other recent works: photometric redshifts for galaxies 
in the Great Observa tories Origins Deep Survey (GOODS; 
[Dickinson et al. 2001) were estimated for 434 galaxies by 
|Mpbashcr ct al. (2004) where in total there were as many as 
18 independent photometric measurements for each galaxy. 
Using the Bayesian method iBeniteall999r they obtained a 
total rms scatter, atot, of 0.11, with an outlier-clipped rms, 
ffred, of 0.047 and an outlier fraction, r/, of roughly 10 per 
cent. This is very similar to the success of the ImpZ code on 
the INT WAS ELAIS sa mple, which had a max imum of only 
5 photometric bands. In lRichards et aP ||200J) . using statis- 
tics based on A^ = {zphot — Zspec), redshifts were estimated 
for 2625 quasars in the SDSS filter system with an overall 
rms, Oz of 0.676 and a^Ed of 0.1 for sources with \A.z\ < 0.3. 

The success of the ImpZ code can also be d i rectly com- 
pared for quasars to the results of iBarger et alj (|2003) who 
got one third of broad-line sources within 25 per cent of 
Zspec . ImpZ fit two thirds of the same CDFN broad-line sam- 
ple within 25 per cent. More recentlv. lKitsionas et al.l (120031) 
estimated photometric redshifts for a sample of X-ray se- 
lected QSOs in the SDSS bands (up to 5 filters), obtaining 
Zphot within \Az\ < 0.3 for 20 of the 30 QSOs at z > 0.4. 
The same statistic can be calculated for INT WAS ELAIS 
'QSO' z > 0.4 sub-sample, where 10 of 20 sources are within 
the same Az tolerance, and for the broad-line and/or com- 
pact CDFN sources at z > 0.4, where 36 of 50 (72 per cent) 
of sources are within \Az\ < 0.3. The ImpZ quasar redshift 
results are clearly highly successful and show that reliable 
redshift estimates for quasars can be achieved. 

It is important to note that the intrinsic variability of 
quasars offsets the magnitudes measured in various bands 
depending on the actual epoch of observation. This is cer- 
tainly an issue for the INT WAS survey which was completed 
over several years. The variability scrambles the real spec- 
trophotometric data taken over an extended period whereas 
the template fits assume that the the photometry in each 
band is taken at the same time (or that there is no vari- 
ability) and is therefore samp ling the same und erlying SED. 
This problem was noted in I Wolf et al.l (|2001j) who found 
around half their quasar sample had completely incorrect 
photometric redshifts, presumably mainly due to this prob- 
lem. Clearly future surveys need to be scheduled with this 
effect in mind. 

The SWIRE survey (in progress) combines both 
ground-based optical photometry and IR data (3.6, 4.5, 5.6, 
8, 24, 70 and 160^m). Virtually all SWIRE galaxies are ex- 
pected to be detected in the IRAC 3.6^.m band. Since little 
contribution from dust is expected at this wavelength, this 
band can be used to improve the photometric redshift es- 
timates. Additionally, at least for low redshift galaxies, an 
estimate can be made of the total stellar mass in the galaxy. 

Where galaxies are detected at 4.5, 5.8 and 8^.m an es- 
timate can also be made of the luminosity in the cirrus emis- 
sion component, after subtraction of the predicted starlight 



contribution. The ratio of hdrr to, say. Lb should be related 
to the A„ value. 

The crucial improvement of a survey such as SWIRE 
over deep fields such as HDF will be the far larger volume 
sampled, so that cosmic variance can be properly addressed. 
Although photometric redshift techniques have been applied 
successfully to the HDF, their true statistical power comes 
from application to large area surveys such as SWIRE, which 
have large numbers of galaxies. Application of the ImpZ code 
to the SWIRE survey will be presented in a future work. 



8 CONCLUSIONS 

In this work, photometric redshifts have been studied us- 
ing SED template-fitting on two spectroscopic redshift cat- 
alogues (INT WAS ELAIS and CDFN). The overall accu- 
racy of the redshift code, ImpZ, is found to be good, with 
Az/{l+z) < 0.1 for 92 per cent of galaxies (where this statis- 
tic is calculated from the 138 'GALAXY' sources in the INT 
WAS ELAIS catalogue). The addition of A„ freedom in the 
fitting improves the redshift solution, suggesting that some 
information about the true A^ of the source is also returned. 
Application to sources with known A^, will help to quantify 
this information (Babbedge 2004; PhD thesis, in prep.). The 
template-fitting method is also extended to quasars via the 
inclusion of AGN templates and is found to be reasonable, 
with Az/{1 + z) < 0.25 for 68 per cent of quasars (based on 
the 25 'QSO' sources in the INT WAS ELAIS catalogue). 
It is noted that the inherent time- variability of quasars re- 
duces the effectiveness of the technique for surveys where 
the photometry is collected over an extended period. 

A combined galaxy-quasar approach to template-fitting 
photometric redshift techniques has been presented and the 
results and analysis presented in this work clearly show that 
photometric redshifts can be calculated with good success 
for both galaxies and quasars, as long as certain consider- 
ations are taken into account: it is important to limit the 
absolute magnitude of the solutions to prevent unphysical 
results; use available optical morphology /stellarity informa- 
tion in order to pre-select sources that are more likely to be 
quasars than galaxies (since we expect quasars to be point- 
like whereas most galaxies will be extended in some way). Al- 
though the application of template-fitting to a single source 
in order to derive its redshift can never be completely relied 
on, due to the inherent degeneracies and limitations of a 
method reliant on perhaps four or five data-points, the sta- 
tistical information drawn from application to a large cata- 
logue of objects is extremely powerful and can be used for 
many different investigations. 

The ImpZ code will be applied to the entire re- 
calibrated INT WAS ELAIS Nl and N2 data in order to 
investigate the evolution of extinction and star formation 
rates in a companion to this paper; Babbedge et al. (2004; 
in prep.). 
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